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Abstract 

Results for da/dfl and in the reaction pp — ► 7r _ 7r + are predicted by a simple quark model. 
They are compared to recent experimental data from LEAR, as well as to previous predictions from 
nucleon-exchange models. At low energy the quark model does better than the nucleon-exchange 
models, but the overall comparison to experiment remains poor. In particular, the double-dip 
structure of the experimental An data is only partly represented. This shortcoming of the simple 
quark model is traced back to a too small J=2 amplitude. This has interesting implications for the 
range of this specific annihilation process. 



PACS numbers: 25.43. +t, 21.30.+y, 13.75.Cs, 12.40.Qq, 12.40.Aa 



1. Introduction 

Large left-right asymmetries in the reaction pp — > 7r~7r + are observed [1-5] when the target 
proton is polarized perpendicularly to the scattering plane. For example, when the spin of the 
target proton is up, the is~ particles prefer to emerge to the right of the p beam-direction at low 
initial antiproton momenta (e.g. plab = 500 MeV/c, or Tlab = 124 MeV), while at high momenta 
(e.g. plab = 1500 MeV/c, or Tlab = 832 MeV) the i\~ particles prefer to emerge to the left. 

The reaction pp — > 7t~tt + is one of the simplest pp annihilation channels because it depends on 
only two helicity amplitudes. In terms of partial waves it depends only on the spin-triplet coupled 
waves, while isospin takes the values 1=0 for even J and 1=1 for odd J. Therefore the existence of 
accurate data, in particular for the analyzing power (also called A y ) is a challenge to microscopic 
models (like the quark model) for the pp annihilation into that particular channel. 

For a detailed review of pp annihilation in general we refer to ref . [6] . Previous work that explicitly 
compares theoretical predictions for the above reaction pp — > 7t~tt + with the available da/dQ and 
An data is based on an annihilation mechanism driven by nucleon- and delta-exchange [5,7,8]. 
Reasonable agreement with the experimental structure of A^ is obtained at antiproton momenta 
in the range of 650-1100 MeV/c (Tlab = 203-575 MeV). The agreement improves towards the 
higher end of this momentum range. However, at lower momenta (plab — 467-585 MeV/c, Tlab — 
110-167 MeV), so far no model has predicted an asymmetry of overall negative sign at all angles, as 
is observed experimentally (see ref. [4,5]). Some reservation exists about nucleon-exchange models. 
Certainly, a legitimate concern of models based on nucleon-exchange is the distant extrapolation of 
the 7rNN form factor, from on-shell nucleons, into a region where one of the nucleons has become 
spacelike. Typically in this extrapolation an effective cutoff- mass of 1.5-1.7 GeV has been used. 

It is therefore of interest to investigate an alternative approach to a microscopic annihilation 
mechanism, which is based on the constituent quark model. In a simple quark model one does not 
specify which object is being exchanged between the anti-proton and the proton, but rather how 
the three quarks and three anti-quarks in the initial state can be rearranged and partly annihilated 
in order to form two qq pairs, representing the final two pions. 

In this paper we study the predictions at low energy for the differential cross section da/df2 and 
the analyzing power A^ if the annihilation mechanism is described by a simple quark model. We 



2 



compare the results of the quark model with two other models where the annihilation mechanism 
is given by nucleon-exchange. From comparing with the experimental data from refs.[4,5], we are 
able to draw some conclusions about the specific annihilation of pp into 7t~tt + . 

In a more general context, one may also foster the hope that because of the short range nature of 
the annihilation process the reaction pp — > 7r~7r + could be a testing ground for possible differences 
between a description in terms of quarks and a description in terms of mesons, nucleons, and isobars. 

2. Annihilation mechanism 

To describe the annihilation of anti-proton proton into two pions we employ here the generalized 
3 P an d 3 5i mechanism of ref. [9] . The proton (and antiproton) wave function is described as a 
Gaussian 

*iv(rL,f 2 ,f 3 ) = ATjv cxp [ - ^ £ (f; - f N ) 2 ] X N (spin, isospin, color), (1) 

where r- x are the quark coordinates and r N is the nucleon coordinate. If we take the pion to be an 
S-wave meson, its wavefunction is given by 

$m (ri,n) = N M exp [ - ^ £ (fi - f M ) 2 X M (spin, isospin, color), (2) 

where r x and r 4 are respectively the quark and antiquark coordinates, and f M is the coordinate of 
the meson. Typical parameter values are a = 2.8 fm -2 and (3 = 3.23 fm~ 2 , giving a nucleon radius 
of 0.60 fm and a meson radius of 0.48 fm. 

The transition potential for this process, predicted by the simple quark model, was worked out 
in ref. [9]. It turns out to be a non-local potential, dependent on the relative pp coordinate r, the 
relative nn coordinate f and the Pauli spin operator a. If one allows for two mechanisms, the 3 Po 
and 3 5'i mechanisms, to be responsible for the annihilation of one of the qq pairs, the potential can 
be written as 

Knn = V( 3 P ) + A V^x). (3) 
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We repeat here the form of V( 3 Pq) derived in ref. [9]: 



\/( 3 P )(r',r-) ~ 
{A v ia.r' sinh(Cr'.f) + By ia.r cosh(Cf ' .r)} exp{A r' 2 + B r 2 ), 



(4) 



where A 



a(a + (3) 



(5) 



v 



2{4a + 3/3)' 



By 



A 



3(5a 2 + 8a(3 + 3(3 2 ) 
2(4a + 3(3) 
a(5a + 4(3) 



(7) 



(6) 



B 



2(4a + 3(3)' 
3{7a 2 + 18a^ + 9(3 2 ) 
8(4a + 3/3) 



(8) 



C 



3a(o; + /3) 
2 (4a + 3(3)' 



(9) 



The expression for V( 3 Si) consists of a longitudinal part and a transversal part. The longi- 
tudinal part has the same general form as eq.(4), the transversal part is obtained from eq.(4) by 
interchanging sinh and cosh. Both expressions, with the appropriate coefficients are given in ref. [9]. 
The parameter A in eq.(3) gives the relative strength of the two mechanisms. The range and the 
nature of the non-localities of V ann is a function of the size-parameters a and (3 of the proton and 
the pion respectively, mentioned in eqs.(l,2). The annihilation potential V ann of eq.(3) contributes 
to all total angular momenta J of the pp system. 

3. Observables and basic amplitudes 

As initial state we choose the pp wavefunction as given by the 1982 Paris pp model [10] including 
the Coulomb interaction. As final state we choose non-interacting pions. Later, we will discuss the 
implications of different cases of initial state interactions and final state interactions, which may be 
substantial [11], in a different paper. 
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The observables, which have been measured in this reaction, are the differential cross section 
da/dfl and the analyzing power An- Both observables are expressed in the two helicity amplitudes 

F ++ (9) and F H (9) which fully describe the reaction pp — > tt~tt + . The angle 9 is the CM angle 

between the outgoing 7r~ and the incoming anti-proton. The observables considered are 

da/cM = ( |F ++ | 2 + |F+_| 2 )/2, (10) 
A N da/dtt = Im (F ++ F*_). (11) 

For completeness and later reference, we note the relation of the two helicity amplitudes to the 
angular momentum amplitudes // and the Legendre polynomials Pj(cos9) and their derivatives 
Pj(cos9). The indices J and £ are respectively the total angular momentum, and orbital angular 
momentum of the pp system. The angular momentum of the 7r~7r + system is £ n7T = J. We have [12] 

F ++ (9) = - vW+lj/2 (Vlfj-! - y/(J+l) fj +1 )Pj(cos9), (12) 
P j 

and 

F + _(9) = I (y/l/Jtf-i + (13) 

p J 

The angular momentum amplitudes // are basically two-dimensional integrals in r' and r of the 
type 

Im = J Jr' 2 dr'r 2 dr<l>y(r')V m (r',r)^%(r). (14) 

Here Vjj^(r', r) originates from the angular momentum decomposition of the transition potential 
Vam^f r, a) and takes on two possible forms 

V lje {r',r) = r exp [A r' 2 + B r 2 } jJ od {Cr'r), (15) 
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and 



V 2M (r',r) = r'ex V [Ar' 2 + Br 2 } jT° d {Cr'r). (16) 

The symbol j™° d m eqs. (15,16) stands for the modified spherical Bessel function of the first kind 
of order £. In eq. (14) above $y(r') and (r) are the final mr and initial pp wavefunctions 
respectively. 

As expected the integral of the above eq.(14) probes the pp wavefunction and 7T7t wavefunction 
at rather short distances since the values of A, B, and C are typically A = -1.80 fm~ 2 , B = -5.59 
fm -2 , and C = -1.21 fm -2 . On the other hand one needs to realize that, for example, the wave 
function ^ p f e (r) vanishes at short distances r. This is caused by the total pp annihilation pre- 
dominantly in channels different from 7T7T, which, in case of the Paris model [10] is described by a 
phenomenological annihilation potential. This of course means to some extent a double counting 
of the specific annihilation into tttt. But since the annihilation into tttt only accounts for about a 
percent of the total annihilation, double counting due to use of the Paris model for the initial pp 
interaction can be ignored. We can safely use this method to study the microscopic behaviour of 
the specific annihilation pp — > tt~tt + . From the above it is clear that integrals of the type of eq.(14) 
are sensitive to the model used for pp initial state interaction, and at the same time sensitive to the 
model for tttt final state interaction. We return to this aspect in a separate paper. 

4. Results 

If we take for ty p / e (r) the Paris model of 1982 [10] and for (r') simply plane waves, (no tttt 
interaction as they go out) , we obtain the result at plab = 497 MeV/c (Tlab = 123 MeV) for 
differential cross section da/dfl and analyzing power An, shown in fig.l for the case A = - 0.5. The 
parameter A which gives the relative strength of the 3 Si versus 3 Pq mechanisms, can be varied. The 
other free parameter, the overall strength off V ann is of no consequence for An while for da/dQ it is 
chosen so as to obtain the correct experimental total cross section. The values for a and (3 are kept 
fixed. With A = -0.5 one is able to obtain an overall negative An at Plab — 497 MeV/c (Tlab — 123 
MeV), although the forward minimum is not deep enough. The relative poorness of the prediction is 
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however still better than the predictions at similar energies by the nucleon-exchange models [5,7,8] 
as is shown in fig. 2. Both N-exchange models give a forward positive value of An at low energy, 
as long as the final tttt interaction is ignored. On the other hand, when tttc FSI is included in 
the N-exchange model [11], the forward An turns negative but the backward An becomes positive 
instead. As mentioned above, experimental data for An at these energies are negative at all angles. 
By itself, it is amazing that this simple quark model is able to predict observables that show any 
resemblance to experiment. 

The number of partial waves needed at this momentum plab is rather small. In fig. 3 it is shown 
that already J<4 and J<5 cannot be distinguished. This means that at this energy J<5 is certainly 
sufficient. 

The prediction at higher momentum, Plab — 679 MeV/c (Tlab — 210 MeV) with no change 
in the parameters A, a, (3, and the overall strength of V ann , is given in fig. 4. Again the forward 
negative dip in An is not deep enough. It is however encouraging that with the same value A = -0.5 
the quark model now predicts a positive value for An near cos# = 0, while the negative forward and 
backward values for An remain. The experimental data at this energy [4,5] show the same trend of 
a double dip with a positive hump at 90 degrees. It is still sufficient at this energy, to include only 
amplitudes with J<5. For easy comparison also the nucleon-exchange predictions of refs. [5,7,8] are 
shown in fig.4. 

One particularly observes, in this low energy range, the general failure of all models to predict 
the large forward dip in An, simultaneously with a large negative value in the backward direction. 

5. Discussion of the results 

We note that predictions at low energy for da/dQ and An in the simple quark model are better 
or of the same quality as for N-exchange models. Of course we have a free parameter A to play with 
that gives the relative strength of the 3 Si and 3 Po mechanisms, as well as an overall factor that sets 
the scale of the cross section. 

At this point one can ask the question, why can one not do better? To answer that question, it 
is helpful to explicitly check the angular behaviour of the terms in the sum over J in the amplitudes 
F ++ and of eqs.(12,13) which are 



7 



F ++ (0) ~ - y/l/2 3 P + y/s/2 ( 3 S 1 - V2 3 D 1 ) cos# 
+ 03/2 ( V2 3 P 2 - VS 3 F 2 ) (3 cos 2 # - l)/2 + . . . , 



(17) 



and 

P+_(#) ~ 0/2 ( 3 5i + 072 3 Pi) sin# (18) 
- 0/2 (y/l/2 3 P 2 + 0L/3 3 F 2 ) (3 sintfcostf) + .... 

From inspection it is clear that for only J=0,1 the asymmetry An has the form 

A N ~ sin(6») (a + ai cos (6»)). (19) 

At low energies, the expression of eq.(19) seems to be the basic form of the predicted A N in both 
N-exchange models, and to some extend the same holds for the quark model (see fig. lb). To get 
an angular behaviour of A N that is of the pronounced double-dip structure as in experiment, one 
needs at least a rather substantial J=2 (e.g. 3 P 2 — 3 F 2 ) contribution. 

This can be clearly demonstrated by a simple toy-model with only 3 P 0? 3 Si, and 3 F 2 pp am- 
plitudes (replacing 3 F 2 by only 3 P 2 does not work as well due to their opposite signs in F ++ (9) 
in eq.(17)). In fig.5 we show the result of such a toy-model where 3 P , 3 Si, and 3 F 2 have relative 
strength 1 : 0.8 : 0.5 respectively. The phases of 3 P and 3 F 2 are the same, 3 S'i has a phase differ- 
ence of 37r/2 with the other two amplitudes. Here the double-dip structure of the asymmetry An, 
shown in fig.5, is due to a strong J=2 presence. 

Going back to the quark model, we observe that indeed there the J=2 amplitude is relatively 
small compared to the J=0,1 amplitudes, and that is the main reason why An has not the pro- 
nounced double-dip structure as seen in experiment. 
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6. Conclusions 

One shortcoming of both quarkmodel and N-exchange models may be that the transition poten- 
tial V ann (r ',f,s) is of too short a range to make the J=2 amplitudes play a significant role at low 
energies. In the simple quark model the range of V ann is determined by the coefficients A,B, and 
C in the exponent. Their values are set by the nucleon size parameter a and pion size parameter 
f3. In order to increase the range of the potential (and enhance the role of the J=2 amplitude) one 
can decrease the values of a and (3. 

As an example we show in fig.6 the result at plab = 497 MeV/c for a = 2.8 fm" 2 and lowering 
(3 to 1.73 fm" 2 (dot-dashed curve). This would correspond with a pion radius of 0.66 fm. The 
previous value was r n = 0.48 fm. Indeed one sees that the double-dip structure of is enhanced. 
At the same time the differential cross section da/dQ is somewhat more forward peaked, in better 
agreement with the experimental shape. In the case of da/dQ the overall strength was adjusted to 
give the same value at cos 6 = -0.5. The corresponding values for A, B, and C are A = -1.79 fm -2 , 
B = - 3.87 fm" 2 , and C = - 1.16 fm" 2 . Similar effects in the asymmetry can be obtained by 
keeping (3 fixed at j3 — 3.23 and reducing a. The dashed curve in fig.6 shows the result for a = 
1.24 fm -2 (r p = 0.90 fm ). In that case the values for A, B, and C become A = - 0.81 fm -2 , B = 
- 4.52 fm" 2 , and C = - 0.57 fm" 2 . For the above two cases the value of A is - 0.45. In both cases 
the role of the J=2 amplitudes is enhanced, when compared with the model with the conventional 
values for a and f3 (solid curve). 

Overall, it seems that experiment requires at this energy a significant J=2 amplitude. One way 
to accomplish this, is to increase the range of the transition potential as shown above. Of course 
one can also search for the reason of the discrepancy between theory and experiment, in the model 
dependence of the initial and final state interaction. Work on this is forthcoming. 

Finally, because the reaction pp — > 7r°7r° has only isospin 1 = amplitudes and therefore only 
even J contributions, the measurement of A^ for that process, if ever performed, may shed addi- 
tional light on the microscopic process of pp annihilation into two mesons. 



9 



Acknowledgment 

One of the authors (W.M.K) wishes to thank the Institut de Physique Nucleaire at Orsay for its 
warm hospitality. The authors are grateful to B. Loiseau for stimulating discussions and for making 
available the pp results of the Paris model. 

References 

[I] E. Eisenhandler et al, Nucl.Phys. B69, 109 (1975). 
[2] A. A. Carter et al, Nucl.Phys. B127, 202 (1977). 
[3] T. Tanimori et al, Phys.Rev.Lett. 55, 1835 (1985). 

[4] R. Birsa et al, Nucl.Phys. (Proc.Suppl.) B8, 141 (1989). 

[5] R. Birsa et al, Proc. 1st Biennial Conf. on Low Energy Antiproton Physics (Stockholm 1990), 
ed P. Carlson, A. Kerek, S. Szlagyi (World Scientific, Singapore, 1991) p. 180. 
[6] C. B. Dover, T. Gutsche, M. Maruyama, A. Faessler, Prog, in Particle and Nuclear Physics 29, 
87 (1992). 

[7] B. Moussallam, Nucl. Phys. A407, 413 (1983); A429, 429 (1984). 

[8] V. Mull, J. Haidenbauer, T. Hippchen, K. Holinde, Phys. Rev. C44, 1337 (1991). 

[9] G. Bathas, W. Kloet, Phys. Lett. B301, 155 (1993). 

[10] J. Cote, M. Lacombe, B. Loiseau, B. Moussallam, R. Vinh Mau, Phys. Rev. Lett. 48, 1319 
(1982). 

[II] V. Mull, K. Holinde, J. Speth, Phys. Lett. B275, 12 (1992). 
[12] A. D. Martin, M. R. Pennington, Nucl. Phys. B169, 216 (1980). 



10 



FIGURE CAPTIONS 



FIG. 1. Differential cross section and analyzing power at plab = 497 MeV/c. Parameter values are 
A = - 0.5, a = 2.80 fm~ 2 , and (5 = 3.23 fm -2 . Data for da/dQ are from ref. [3], and A N data are 
from ref. [5] . 

FIG. 2. Analyzing power at plab = 497 MeV/c. The solid curve is the quark model. The nucleon 
exchange model is represented by the dashed line ref. [5,7], and the dotdash line from ref. (8). All 
models are without nn FSI. Data are from ref. [5]. 

FIG. 3. Effect of various partial waves at Plab — 497 MeV/c. The solid curve is the prediction for 
J<5, and cannot be distinguished from J<4. the dashed curve is J<3, and the dotdash curve is J<2. 

FIG. 4. Differential cross section and analyzing power at plab = 679 MeV/c. The solid curve is the 
quark model, with the same parameters as in fig.l. The dashed curve is ref. [5,7], and the dotdash 
curve is ref. [8]. Data for da/dQ are from ref. [3], and A N data are from ref. [5]. 

FIG. 5. Result from simple toy-model with 3 P , 3 S±, and 3 F 2 amplitudes. For details see text. 

FIG. 6. Results at p LA B = 497 MeV/c for a = 2.8 fm~ 2 and (3 = 1.73 fm~ 2 (dotdash curve). Dashed 
curve is for a = 1.24 fm -2 and (3 = 3.23 fm~ 2 . In these first two curves A = -0.45. The solid curve 
serves as reference to the earlier values a = 2.8 fm -2 and (3 = 3.23 fm -2 , where A = -0.50. Data 
are from ref s. [3,5]. 
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